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PLASTIC BUCKLING OF A RECTANGULAR PLATE UNDER EDGE THRUSTS 

By G. H. Handelman and W. Prager 


SUMMARY 

The fundamental equations for the plastic buckling of a 
rectangular plate under edge thrusts are developed on the basis 
of a new set of stress-strain relations for the behavior of a metal 
in the plastic range. These relations are derived for buckling 
from a state of uniform compression. The fundamental equa- 
tion for the buckling of a simply compressed plate together with 
typical boundary conditions is then developed and the results 
are applied to calculating the buckling loads of a thin strip, a 
simply supported plate, and a cruciform section. Comparisons 
with the theories of Timoshenko and Ilyushin are made. Fin- 
ally, an energy method is given which can be used for finding 
approximate values of the critical load. 

INTRODUCTION 

This paper is concerned with the plastic buckling of a 
rectangular plate which, previous to buckling, is under a 
uniform compressive stress tr 0 in the direction of one of its 
edges. In the case of elastic buckling in which <x 0 remains 
below the clastic limit of the plate material, it is well known 
that the buckling stress depends on the dimensions of the plate 
and on the manner in which it is supported (cf . reference 1 , ch. 
7). In the case of the plastic buckling of beams, on the other 
hand, Engesser (reference 2) and Von Karman (reference 3) 
developed a satisfactory theory based on the fact that for a 
fiber which is compressed beyond the elastic limit the tangent 
modulus (i. e., the ratio of the variation of strain to the cor- 
responding variation of stress) assumes different values depend- 
ing on whether the variation of stress constitutes an increase 
or a relief of the existing compressive stress. 

Generalization of this theory to the pla'stic buckling of 
plates has repeatedly been attempted. These attempts can 
be divided into two groups which may be labeled formal and 
analytical generalizations. The formal generalizations start 
from the remark that the formulas of the Engesser-Von 
Kerman theory of the plastic buckling of beams differ from 
the well-known formulas for the elastic buckling of beams 
only by the fact that the so-called “reduced modulus” 
replaces Young’s modulus. A formal generalization of the 
Engesser-Von Karman theory to the plastic budding of plates 
is therefore obtained by introducing the reduced modulus 
into the formulas for the elastic buckling of plates in such a 
manner that the results of the Engesser-Von Karman theory 
are obtained in the case of a narrow rectangular strip which 
is free on its long edges and simply supported on the short 
edges where it carries a compressive load. Of course, this 
formal generalization is more or less arbitrary and leads by 


no means to a unique result. Formulas of this type have 
been suggested by Bleich (reference 4, p. 216 ff.) and Timo- 
shenko (reference 1, p. 384). 

In contrast with these formal generalizations of the Enges- 
ser-Von K&rm&n theory, the analytical generalizations do not 
merely introduce the reduced modulus of the theory of beams 
into the formulas for the elastic buckling of plates. Instead, 
the analytical generalizations go back to the considerations 
by which the reduced modulus is derived and try to apply 
these to the case of a buckled plate. Generalizations of 
this kind have been previously presented by Kaufmann 
(reference 5) and Ilyushin (reference 6). As is shown in the 
present report, however, these authors use stress-strain 
relations which do not fulfill certain postulates of the theory 
of plasticity; the correctness of their results must therefore 
be. questioned. 

The present paper aims at developing a theory of the 
plastic buckling of plates which takes full account of the 
modern theory of plasticity. The stress-strain relations in 
the plastic range are discussed at considerable length in the 
first section of the ANALYSIS, and it is shown that, for an 
adequate treatment of buckling phenomena, a theory of 
plastic flow is indicated rather than a theory of plastic defor- 
mation of the type used by Kaufmann and Ilyushin. The 
precise definitions of these terms and the basic considerations 
suggesting the use of a theory of plastic deformation for 
problems such as buckling are fully discussed in the ANALY- 
SIS. A particular theory of plastic flow suitable for the treat- 
ment of the problems under consideration is developed in the 
first section and its relations with other theories of plasticity 
are pointed out. It is shown that in the particular case of 
a plate buckling out of a state of simple compression there is 
very little freedom in the choice of the stress-strain relation 
if it is to fulfill certain simple postulates. This means that 
all the empirical information which is necessary for the theo- 
retical treatment of the plastic budding of a rectangular 
plate under edge thrusts can be obtained by a simple com- 
pression test. 

The second section presents the development of the funda- 
mental equation of the plastic budding of a simply compressed 
plate, and the appropriate equations describing typical bound- 
ary conditions are given in the third section. The remaining 
parts contain several examples, which are carried out in 
detail, as well as an equivalent energy principle which proves 
to be very useful for approximate computations. Finally, 
the appendixes contain detailed discussions of several tech- 
nical points raised in earlier parts of the paper. 
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SYMBOLS 
length of plate 

constant in compressive stress-strain law 
expression in variational principle 

jcoefficients in plastic stress-strain law 

width of plate 

expression in variational principle 


C= (X — 1) /[ (5 — 4r)X — (1 — 2r) 2 ] 


Cl, C 2 , Cz, Ci 


D 

Dll, D\2, D 2 2 


arbitrary constants appearing in equa- 
tion for w 

flexural rigidity of plate h 3 E u / (1—p 1 )^ 

coefficients in plate equation for plastic 
flow 


Dn'=Dn/D,D 12 '=D l2 /D,Dv f =D 2i /D 

Dn*,Di 2 *,D 22 * coefficients in plate equation for Ilyu- 
shin’s theory of plastic deformation 


Du=Dn*/D, Di 2 =D t2 *jD, D 22 =D 22 *ID 

E tangent modulus in compression 

E 0 Young’s modulus 

E* Von Karman’s reduced modulus 

E s secant modulus obtained from compres- 

sive stress-strain diagram 
section of buckled middle surface for x= 
Constant 
thickness of plate 
moment of inertia of cross section 

cr 0 a 2 (l — v 2 ) 
ir 2 h 2 E 0 


fiy) 

h 

I 

k—12 


K= (2 — v)Ki J r{2v — l) K 2 

K\ = d 2 vj/dx 2 ,K 2 = b 2 w/by 2 ,K i2 = 2d 2 ib/dxdy 


Nx 

N v 


rate of change of stress resultant in x- 
direction 

rate of change of stress resultant in y- 
direction 


P 2 = 


( m7r \ Dn 
v a J D 22 


P total compressive force ( <r 0 bh ) 

Q=[D \ 2 — (1 — p)D\/D 22 


2 _ooh_/ mi r\ 2 D\\D 22 —D\2 
^ D \ 2 \ a ) Di 2 D 22 

r — VjP(2+P) 

R side ratio ( b/a ) 


s= V?( 2 — P) 

t time 


u= 8r 2 + 12^—23 

w deflection rate 

x, y, z rectangular Cartesian coordinates; x,y- 

plane coincides with middle surface of 
unbuckled plate 


2o 


(2 — r)ei+(2r— l)t 2 

k 


a ^ c( 5 — 4p) 


a 


& 


ratio of Von Karman’s modulus to Young’s 
modulus (in section “Buckling of a 
simply supported plate” only) 

constant in compressive stress-strain law 


8 = 


A 


2 ^f 2 ^ 0+ +2 

D 11 D 22 ' — D X2 ' 2 

d 22 2 


for K> 0 
for 0 


m 

M x 

My 

Akxy 

n 

N x 


number of half waves in buckled configur- 
ation 

rate of change of bending moment about 
(/-axis 

rate of change of bending moment about 
2 -axis 

rate of change of twisting moment 
integer 

reduced compressive stress resultant 

(a 0 h/E 0 ) 


€ 


uniaxial strain 


dc x j djty, dz 2, 

dy X y, dy UZj dy zx 

dt x , dty , 
de z , dy X y 


) infinitesimal strain increments present in 
) buckling 

| reversible (elastic) strain increments 


de x ", dty", 
de z ", dy xv " 


| permanent (plastic) strain increments 



PLASTIC BUCKLING OF A RECTANGULAR PLATE UNDER EDGE THRUSTS 


3 


«i(a, y) 


normal strain rate in middle surface in in- 
direction 


h(x,y) 

fo=22oA 

fo + = a+ V« 2 “ ’I 

fo~ = -a-V« 2 -l 

rf-sb 


normal strain rate in middle surface in y- 
direction 

shear strain rate in middle surface 


v 

%=rb 


ratio of Von Karmfm’s modulus to 
Young’s modulus ( E*/E 0 ) 

ratio of Young’s modulus to tangent 
modulus ( E 0 /E ) 

Poisson’s ratio 


uniaxial stress 


or 

Oct 

<r 0 

&Xj Gy) G z 
T x\ f) T yz) T zx 
d & X) d 0 ” y j d (T 2 ) 

dr xy , d-T yz) dT 2X 

'!'—(*> ( l — 9 


CO 


co' T=dw/do t 
D = co/(l — co) 


• dQ, . 

12= j <7* 

CtC i 

O' =dt2/dcr i 


+ 


intensity of stress (Voy'+ay’— o x c v + 3 t x ~ 2 ) 
critical compressive stress 
original compressive stress in plate ( — a x ) 
normal stress components 
shear stress components 
) infinitesimal stress increments present in 
] buckling 



function of intensity of stress cr. 


superscript denoting values on unloading 
side of neutral surface 
superscript denoting values on loading 
side of neutral surface 


ANALYSIS 


STRESS-STRAIN RELATIONS FOR BUCKLING FROM A STATE OF 
UNIFORM COMPRESSION 


The mechanism of buckling beyond the elastic limit is 
relatively complicated because the material, which was 
originally in a state of simple compression, is loaded in some 
regions and unloaded in others during the buckling process. 
Consequently, the stress-strain relations must be considered 


in some detail with special reference to the problem of load- 
ing beyond the elastic limit followed by unloading. 

The material must exhibit strain-hardening if the deter- 
mination of the buckling stress is to constitute a problem. 
Indeed, for a perfectly plastic material which yields under 
constant stress, Von K&rmdn’s reduced modulus vanishes 
once the initial compressive stress has reached the yield 
limit. This means that the bending stiffness is reduced to 
zero and buckling must be expected quite independent of 
the dimensions of the bar. 

Stress-strain laws for materials which exhibit strain- 
hardening can be divided into two types which, for con- 
venience, will be called “theories of plastic deformation” 
and “theories of plastic flow.” According to the first group, 
there exists a one-to-one correspondence between stress and 
strain in the plastic range, as well as the elastic, provided 
that the material is being loaded. The stress-strain law of 
the well-known Hencky-N&dai theory (reference 7, ch. 14, 
and reference 8) and the law used by Ilyushin (reference 6) 
in his discussion of plastic buckling are typical theories of 
plastic deformation. On the other hand, the theories of 
plastic flow are based on the assumption that, for a given 
state of stress, there exists a one-to-one correspondence 
between the rates of change of stress and strain in such a 
manner that the resulting relation between stress and strain 
cannot be integrated so as to yield a relation between 
stress and strain alone. Typical examples of theories 
of plastic flow are the stress-strain relations developed 
by Prager (reference 9) and Handelman, Lin, and Prager 
(reference 10). A particularly important difference between 
these two basic theories of plasticity lies in the fact that the 
strain which corresponds to a certain state of stress, accord- 
ing to the theory of plastic deformation, is entirely inde- 
pendent of the manner in which this state of stress has been 
reached, whereas, according to the theory of plastic flow, the 
strain depends on the manner in which the state of stress is 
built up. 

The stress-strain relations to be used in the analysis of 
the plastic buckling of a rectangular plate under edge thrusts 
form a special case of those developed by Handelman, Lin, 
and Prager in reference 10. In this particular case, how- 
ever, it is possible to develop the stress-strain relation in a 
quite elementary manner, and the inherent difficulties of 
the theories of plastic deformation can be seen from a 
slightly different point of view. It appears worth while, 
then, to examine these relations in some detail with special 
reference to the problem which forms the subject of the 
present report. 

In the following, the stresses and strains in the budded 
plate will be referred to a fixed system of rectangular Car- 
tesian coordinates x, y, and z. The r,?/-piane of this coor- 
dinate system coincides with the middle surface of the un- 
buckled plate, and the axes of x and y coincide with two of 
its edges, the other edges falling on the lines x=a and y—b. 
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Figure 1.— Plate under uniform compressive stress in direction of z-axis prior to buckling. 


Prior to buckling, the plate is under a uniform compressive 
stress <r 0 in the direction of the z-axis (fig. 1). The investi- 
gation of the stability of the state of stress 



requires the knowledge of the relations between the infinites- 
imal increments of stress do x , do v , do z , d.T xy , dr V2 , and dr zz 
and the corresponding increments of strain dt x , dt y , dt z , 
dy xy> dy vz , and dy 2X . Within the framework of plate theory, 
a z =T xz =T yz = 0, even in the buckled state, and hence do x ~ 
drxz—drj 2 =0. Accordingly, dy xz =dy V2 =0. Within the elas- 
tic range the remaining increments of stress and strain are 
related to each other by means of 

E a de x =do z — vda v 

E 0 de v — vdo z I da y 

y ( 2 ) 

E 0 de x = — vdd x — vdi t v 
E yly x ,j~ —0 I d)d,i' xv ^ 

where E a denotes Young’s modulus and v, Poisson’s ratio. 
Before an analysis of the plastic buckling of the plate can 
be attempted, the relations replacing equations (2) in the 
plastic range must be known. In order to establish these 
relations, it will be convenient to think of the strain incre- 
ments^as consisting of reversible (elastic) and permanent 
(plastic) components: 

dt x =dc x r drdtx" 

dty = dty' -\-dty ' 

}- (3) 

dt z — dt z — j dt 2 

dyxy~dyxy ~\~dy xy ^ 

Primes and double primes denote elastic and plastic com- 
ponents, respectively. The elastic increments of strain are 
related to the increments of stress by means of equations 
(2), in which the left-hand sides must all be written with 
primes now: 


E 0 dt x — do ■ vdoy 

EyfLty vdo x | d(Ty 

r ( 4 ) 

E a de z r = —vdo x — vdo y 

E 0 dy xv ' = 2(1 + v)dr x y ^ 

These relations may be regarded as the definitions of the 
elastic increments of strain. The purpose of the following 
discussion is to establish similar relations for the plastic 
increments of strain. 

The elastic increments of strain, equations (4), depend 
only on the increments of stress and are independent of the 
existing stress o 0 . Moreover, a reversal of the signs of all 
increments of stress leads to a mere reversal of the signs of 
all elastic increments of strain. The plastic increments of 
strain, however, do not have these properties; since they 
must vanish as long as <r 0 remains below the elastic limit, 
they cannot be independent of the existing stress o 0 . More- 
over, if for a given value of o 0 certain stress increments pro- 
duce plastic increments of strain, stress increments of the 
same magnitudes but opposite signs do not produce any 
plastic deformation. In other terms, beyond the limit of 
elasticity an infinitesimal change of stress may be classified 
as loading the material or not according to whether it is 
accompanied by permanent deformation. Infinitesimal 
changes of stress which do not load the material may be 
classified in turn as unloading or neutral. Unloading brings 
the material into a state of stress such that all sufficiently 
small further changes of stress are accompanied by elastic 
deformations only. These basic differences in loading and 
unloading appear somewhat more clearly if the simple 
example of a uniaxial state of stress and strain (say a tensile 
test) is considered. Let o denote the stress and e the strain 
in figure 2 and suppose the material is loaded to the point 
P. The stress-strain diagram for unloading is a straight 



Figure 2.— Stress-strain diagram for loading and unloading for uniaxial state of stress and 

strain. 
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line PA with the same slope as the loading curve at the 
origin 0. The permanent strain corresponding to loading 
up to the point P is measured by OA. Suppose now that, 
after the point P has been reached, the test specimen is 
further loaded to the point Pi ; by this, the permanent strain 
is increased by the amount AA\. In other words, the change 
from P to Pi constitutes loading in the sense just defined. 
On the other hand, if the new state of stress and strain is 
given by the point ft, that is, if the stress has been reduced 
below that at P, the permanent strain is left unchanged. 
Furthermore, any small change of stress from the point P 2 
(strictly, all changes within the ranges P 2 P and Pa A) pro- 
duces an additional deformation which is purely elastic.. 
The material has thus been unloaded. 

For uniaxial stress any change of stress constitutes either 
loading or unloading. A third possibility, designated as 
“neutral” change of stress, exists in the case of combined 
stress. A neutral change of stress, while not accompanied 
by a permanent deformation, brings the material into a 
state such that there exist certain further changes of stress 
which are arbitrarily small and yet produce a permanent de- 
formation. This third condition is illustrated in the analysis 
of the buckling of a plate. It is precisely the possibility of 
the occurrence of neutral changes of stress which distinguishes 
the present problem from that treated by Engesser and Yon 
Karman, for in their case the stress is uniaxial in the budded 
state as well as in the unbuckled. Accordingly, a change of 
stress can bo only an increase of the existing compressive 
stress (loading) or a decrease (unloading). The situation is 
more complicated in the case of a plate. 

Since there is no permanent deformation accompanying 
neutral changes of stress or unloading, de x " =de v ' ' =dt z " = 
dy xv "= 0; and the relations of equations (4) define the total 
change of strain. For loading, however, equations (4) must 
be supplemented by equations of the form 

E 0 dt x " =a 'd<r x +b'da u 


E 0 de v " =a"dcr x - J rb"d<r v 
E 0 de z " =a'''d<j x +b"'dc v 


(5) 


side. For instance, the coefficients of d<r x in the second and 
third equations are equal, as are the coefficients of da v in 
the first and d<r x in the second equations. Which of these 
symmetries, if any, will be maintained in equations (5)? 
The existing state of stress singles out the z-axis, but it 
does not matter which of the other two axes is labeled y 
and which z. Accordingly, 

(9) 

In view of equations (7) and (9), 

a " = a’"----a' ( 10 ) 

These coefficients can easily be expressed in terms of the 
so-called “tangent modulus” corresponding to the compres- 
sive stress tr 0 . Application of equations (4) and (5) to simple 
compression in the ^-direction yields (with da v = 0) 

E 0 dt x =E 0 {de x +de x ") 


= d CT x "j - CL r (1 0 x 

(ID 

or 


jyt d(?x ^0 

de x 1 +a r 

(12) 

where E denotes the tangent modulus. With 


\=E 0 /E 

(13) 

equation (11) gives 


a/ = X— 1 

(14) 

Hence, according to equation (10), 



(15) 


Next, the criterion for neutral changes of stress must be con- 
sidered. Any given infinitesimal change of strain can be 
decomposed in the following manner: 

de x —-g (dtxAdiy-pdtz) +tj (2 dt x dt y de z ) J 


E 0 dy xv " = 2c'dr xv J 

where the coefficients a', b', a”, b", a'" , b"' , and c' depend 
on the existing stress <t 0 . 

As is customary in the theory of plasticity, the plastic de- 
formations will be supposed to represent a mere change in 
shape but no change in volume. Accordingly, 


dt y — o' (di x -\-dtyA'de^) -fx ( ds x A 2 A,, de z ) 


de z =g- (de x -\-de y -\-de^) Tg - ( — de x de v -\- 2 de z ) j (16) 
d.y xv =0-\-dy XV 


de x "+dey" + de/' = 0 


( 6 ) 


dy vz =Q-\-dyyz 


T his relation must hold independently of the values of da x 
and d(T v . Thus, 


a'-f a ,/ +a" , =0 

(7) 

b'+b"+b"'=0 

(8) 


The elastic formulas, equations (4), exhibit a certain 
symmetry of the coefficients appearing on the right-hand 


dy zx =Q) J rd'y zx ) 

The change of strain defined by the first members of the 
right-hand sides of these equations is a uniform expansion 
(or contraction) in all directions. Such a uniform expan- 
sion changes the volume but not the shape of the element 
to which it is applied. The change of strain defined by the 
second members of the right sides of equations (16), on the 
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.other hand, affects the shape of the element but preserves 
its volume. The work done by the existing stress or stresses 
on the change of strain, equations (16), consists of the work 
done on the change of volume represented by the first mem- 
bers of the right sides of equations (16) and the work done 
on the change of shape represented by the second members. 
Since all changes of volume are supposed to be of an elastic 
nature, it seems natural to speak of loading or unloading 
according to whether the work dW which the existing stresses 
do on the change of shape alone is positive or negative. 
Vanishing of this work must then be interpreted as indicating 
a neutral change. 

In the case under discussion, the only existing stress is 
<r x = — (T 0 , and the criterion for loading or unloading is fur- 
nished by the sign of 

dW=--. 3 &de x -dt v -de z ) (17) 

while neutral changes are characterized by 

2 df x — de v — de z = 0 (18) 

Now, for unloading, the entire change of strain is of an 
elastic nature and equations (2) apply. Equation (17) is 
therefore equivalent to 

E a dW= — 7 r [2 (da x — v da v ) — ( — v da x -\-da v ) — ( — v da z — v da v ) ] 

(19) 

Since <r 0 j>0, this expression will be negative, whenever 

2da x —da v ^>0 (20) 

This inequality, equation (20), is thus seen to constitute the 
criterion for unloading. Similarly, the criterion for neutral 
changes of stress is found to be 

2da x — da y =0 (21) 

Changes of stress which satisfy neither equation ( 20 ) nor 
( 21 ), that is, changes of stress for which 

2da x —da v <^0 (22) 

must therefore constitute loading. Another definition for 
the criterion for the three types of change of stress, which 
is found by combining equations ( 20 ), ( 21 ), and ( 22 ), is that 
the change of stress is classified by the sign of the increment 
in the second invariant of the stress deviator, which measures 
the intensity of stress. A detailed account of this alterna- 
tive formulation is found in reference 10 . 

By a suitable choice of the values of da x and da v , the 
expression 2 da x —da v can be made to fulfill the following 
inequalities: 

0>2 da x —da v y — e (23) 

where e is an arbitrarily prescribed small positive number. 
All changes of stress satisfying equation (23) constitute 
loading and are therefore accompanied by plastic deforma- 
tions in accordance with equations (5). For e-> 0 , however, 


these changes tend toward neutral changes of stress for 
which there are no plastic deformations. Furthermore, there 
are no plastic increments of strain when 2da x —da v ^>0. It 
is to be expected that the total strain increments will be 
continuous in the region which marks the transition from 
unloading through the neutral state to loading. Although 
such a statement does not follow specifically from the equa- 
tions of equilibrium or compatibility, continuity should be 
expected in the strain increments. With this assumption, 
the plastic increments of strain, equations (5), should vanish 
whenever the increments of stress satisfy equation (21). 
This furnishes the conditions 

a' + 2b' = 0 
a"+2b"=0 

r (24) 

a'" + 2b'"=0 
c'=0 

Together with equations (14) and (15)-, these equations de- 
termine all coefficients appearing in equations (5), which 
therefore take the form 


E 0 dt x "=(\—l)d 

Eode,” — ^ 

E 0 de z "=-^ 
E 0 dy xu " = 0 


X~l A 

dO'y 


d(T x~\~ A d(T y 


da x 


4 

X-l 


4 < ^ a ' v 


(25) 


It is interesting to note that here again the coefficients of 
da x in the second and third equations are equal, as are the 
coefficients of da y in the first and of da x in the second equa- 
tion. Whereas in the elastic case this type of symmetry in 
the stress-strain relations is a consequence of the isotropy 
of the material, this is no longer so in the case of equations 
(25). Indeed, the equality of a" and a'” (see equation 
(15)) follows from the assumption that the plastic deforma- 
tions do not involve a change in volume. The equality of 
b' and a", on the other hand, might be described as almost 
accidental, the value of the ratio a" la 1 being fixed by the 
assumption just mentioned, while the value of the ratio 
b' la' is fixed by the form of the condition for neutral change 
of stress. 

Combination of equations (4) and (25) finally yields the 
stress-strain relations which will be used throughout this 
paper: 


Eftdt x — \da x ^r-{- 2 ^ da y 

E 0 d^-(v+~-') d<Jy | 

E 0 de z =— + da x — (v — ^ ^ ^ day 


E 0 dy xv =2(1 v) dr x 


( 26 ) 
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It seems worth while to stress once again the assumptions 
on the basis of which these stress-strain relations are derived. 
These are 

(1) Plastic deformations do not involve a change of volume 

(2) The criterion for loading or unloading is furnished by 
the sign of the work dW which the existing stresses do on the 
change of shape produced by the increments of stress 
The first assumption is commonly made in the theory of 
plasticity (cf. reference 7, p. 10) and is confirmed by the 
experiments of Bridgman (reference 11, p. 166). The sec- 
ond assumption is a slight generalization of a similar as- 
sumption which Prager (reference 10) introduced in the case 
of incompressible plastic materials ; more recently, it has been 
used by Ilyushin preference 6) . 

It is interesting to note how far the stress-strain relations, 
equations (26) , differ from those used in previous work on the 
plastic buckling of plates. In the present notation, Kauf- 
mann’s stress-strain relations (reference 5) are 


E$d 6 x — XrZtTj vd(Ty 
E Q de„= — 'Kvda x -\-d<x u y 
E 0 dy X y— (1 + A) (1 -f -v)dr xv „ 


(27) 


(The expression for de z is not given because this strain com- 
ponent is not necessary for the determination of the bending 
and twisting moments in the buckled plate.) It is seen that 
here the coefficient of da v in the first equation and that of 
da x in the second arc unequal. In an earlier paper on the 
plastic buckling of cylindrical shells (reference 12, footnote 
1, p. 422) in which similar stress-strain relations were used, 
Kaufmann comments on this lack of symmetry, recommend- 
ing that the stress-strain relations, equations (25), be checked 
by experiment. Since this typo of symmetry in the present 
stress-strain relations, equations (26) , has been characterized 
as almost accidental, the lack of symmetry in Kaufmann’s 
relations hardly constitutes a sufficient reason for discarding 
the stress-strain relations, equations (27). It is not difficult, 
however, to show that these relations correspond to an un- 
acceptable condition for neutral changes of stress. Indeed, 
subtraction of the elastic increments of strain, equations (4), 
from the total increments of strain, equations (27), yields 
the following plastic increments of strain: 


E 0 de x " = (A — l)da x 

E a dty" = v (\ — 1 )dc x 

E 0 dy xv " = (X— 1) (1 +v)dT X y 
These plastic increments of strain vanish if 


dc r x =0 
dr xv — 0 


(28) 


(29) 


According to Kaufmann’s stress-strain relation, neutral 
changes of stress are characterized by the two conditions 
given as equations (29). If the most general change of stress 
considered here is represented by a point with the coordinates 


d<r x , da-y, and dr xy in a three-dimensional space, the condition 
of equation (21) represents a plane through the origin which 
separates the “region of loading” from the “region of unload- 
ing.” Equations (29), however, define a straight line which 
does not mark off two such regions. 

Ilyushin (reference 6) considers an incompressible material 
and assumes the stress-strain relations for loading to have 
the form 

2<Tz — a - y 


Eot x - 


2(1 — co) 


is _ 2 <Ty <r x 

/t°«„- 2 (i_ a) ) 

3t X y 


E^y X y 


1 — CO 


(30) 


where co is a function of the intensity of stress or defined by 

<r«= Vov + oy*— o x Gy-\-3r X y (31) 

For loading, the increments of stress and strain are then 
connected by 


E a dt x — ^ ( 2 da x da v 


( 2 


E { 


2(1— co) 

'°dq,=2(jrr^ (2 d<? v — 

E 0 dy xu =y^2 1 (dr X y— du) 




d ffx - 2 -^— x dco) 

1 — 07 / 


(32) 


For budding from a state of uniform compression <r x = — a Q , 
in particular, 

dco=^do ( =-y (2da x -da v ) (33) 


where 


co' =dco/dcTj 

Equations (32) then reduce to 


(34) 


E 0 de x — . 


1 


(1 — co) 


- 2 (1 — co — co’o-q) (2da x — da u ) 


E^dcy - 


2(1— co) 
3 

E(jdy x y = X — dr xy 

I — CO 


1 2 | ( 1 — co + 2a 0 ai')da x -j- 2(1 — co) — 11 


-J da v 


1 


(35) 


For unloading, the relations, equations (2), are supposed to 
hold with r=l/2 on account of the assumed incompressibility 
of the plate material: 


Eodtx — g (2 da x d<Xy) 

i?ode„=- 2 - (2(7 (j y da z ) 

E()dy x y == 3dT X y 


(36) 


As to the criterion for loading and unloading, this is again 
supposed to be given by the sign of the expression 


d W — a x d^ x {— a y dt y — | t x ydy x . 


(37) 


W.0621 — 50 2 
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In particular, it is given by the sign of 

— Vo de x (38) 

in the case of budding from a state of uniform compression 
<r x = — <r 0 - Id view of the first of equations (35), this means 
that neutral changes of stress are again characterized by 
equation (21). It is easily seen, however, that for 2 d<r x — 
da v = 0, equations (35) and (36) do not give the same incre- 
ments of strain. Ilyushin’s stress-strain relations are thus 
seen to exhibit an objectionable discontinuity along the 
surface 2d<r x —d<r v =0 which separates the region of loading 
from the region of unloading. A more detailed analysis of 
the effect of this discontinuity in the case of a budding plate 
is found in appendix A. 

FUNDAMENTAL EQUATION OF PLASTIC BUCKLING OF A 
SIMPLY COMPRESSED PLATE 

The technique used in the derivation of the fundamental 
equation of the plastic buckling of a simply compressed plate 
is quite similar to that needed for the same problem in the 
elastic range. (See, for example, the more general problem 
of combined bending and compression of elastic plates in 
reference 1, p. 302.) There is one essential difference, how- 
ever, in that the stress-strain relations given in equations (26) 
must be used in the regions of loading rather than generalized 
Hooke’s law. Consequently, the middle plane of the unbent 
plate will no longer play the role of the neutral surface in the 
budded position. Once the position of the neutral surface 
has been found and the bending and twisting moments deter- 
mined as functions of the second derivatives of the deflection 
of the plate, the equilibrium conditions and the final differen- 
tial equation can be derived in exactly the same fashion as 
that used by Timoshenko in reference 1. 

It will be found more convenient, in the following discus- 
sion, to use “reduced stresses” rather than actual stresses, 
that is, stresses reduced by dividing the actual stress by 
Young’s modulus E 0 . No new notation will be employed to 
denote these reduced stresses ; therefore, care must be taken 
in interpreting the results obtained here in terms of the 
known facts for elastic budding. An attempt will be made 
at such points to keep the notation clear. In addition, the 
use of differentials of stress and strain may lead to some con- 
fusion in deriving the equations of equilibrium for an element. 
Since the stress-strain relations given in equations (4) and 
(26) are linear in these differentials, both sides of the equa- 
tions may be divided by d<> 0, where t may be regarded as 
the time. It should be noted that t appears homogeneously; 
that is, the time scale may be arbitrarily distorted without 
changing the equations. If differentiation with respect to t 
is denoted by a dot, equations (4) and (26) can be rewritten 
as reduced stress-strain relations, for W ^ 0, 



and for W^O, 

i x =ix~~vb v 

Cy V(J x 1 (7y 

i — V(7 x — V(T y 

7X1/ = 2(1 T P) T X y 


(40) 


where W=dw/dt. 

The stress rates b x and b v can be found in terms of the 
corresponding strain rates for loading by solving the first two 
of equations (39). Thus, 


° x ~ (5 — 4r)A— (1 — 2p) 2 t ( x + 3) <k+ 2 (X s - 1 + 2p) i v \ 
<r " = (5'— 4 p) X— (1 — 2v) 2 1 2 ( X “ 1 + 2v ) ^+ 4X «p] 


V (41) 


The criterion for loading 


2(7 x (J y 0 


(42) 


can then be written as 


(5— 4p)X— (1 — 2p) 2 ^ 2 v ) e x+(2v IKRO (43) 

Now, Poisson’s ratio p satisfies the inequality— 1 1/2 (cf. 

reference 13, p. 104); in addition, XS; 1. Consequently, the 
expression appearing outside the brackets is always positive 
and the inequality, equation (43), can be replaced by 

(2-p)« I +(2p-l)ey<0 (44) 


The strain rates appearing in equations (39) and (40) must 
now be evaluated. The strain rates in the middle surface 
will be denoted by ii—ii (x, y), the normal strain rate in the 
z-direction ; e 2 =e 2 (x, y), the normal strain rate in the y- 
direction; and y=y(x, y), the rate of shear strain. Points on 
the normal to the undeformed middle surface arc assumed to 
remain on the normal of the bent middle surface. This 
implies that the strain rates e x , e y , and y xv at any point of 
the plate can be written in the following form: 


ix^ii — zKx ' 
ey — i 2 — zK2 
ixv—y—zKu* 


(45) 


The quantities K u K 2 , and K 12 appearing in equations (45) 
are defined in terms of the rate of deflection w=w{x, y) of the 
middle surface in the following way: 




b 2 w 

b 2 w 

by 2 




2 — 2 


b 2 w 


(46) 


y XV ~~ 2(1 — 1 p) 7 X y 


bxby J 


PLASTIC BUCKLING OF A RECTANGULAR PLATE UNDER EDGE THRUSTS 


9 


Geometrically, and K 2 represent the rates of curvature 
of the middle surface in the x- and y-directions, respectively, 
whereas K n represents the rate of relative twist. The cri- 
terion for loading, equation (44), can now be rewritten in 
terms of the strains of the middle surface and the quantities 
Ki, K 2 , and K u . It is seen that loading takes place provided 

■ (2-v)e l +(2v-l)e i <tz[(.2-v)& 1 +{2v-l)&] (47) 

With the word “sign” to denote the sign of the quantity within 
the parenthesis and vertical bars to denote the absolute 
value of the enclosed expression, the last inequality may be 
transformed into 


z sign [(2— v)Ki+(2v— 1)K 2 ]> 


(2 — v ) (2r — 1) 

I (2 — v)K\-\- (2v — -V}K 2 \ 


(48) 


This inequality can be simplified further by introducing two 
new quantities K and z 0 defined by 


K= (2 - v)K, + (2r - 1 )K 2 ' 

.. {2-v)i l +{2v-\)i 2 y 

2o “ £ 

y 

The inequality, equation (48), becomes then 
z sign (K)>Kz 0 /\]£\=Zo sign ( K ) 


(49) 


(50) 



(52) 


As indicated in equations (51), two cases must be considered 
according to whether £T> 0 or K<^0. For Ky> 0, direct 
computation shows that 


N x — v N y = a +ck(^~ z 0 J 
N v —vN x =i 2 h -i cK - z 0 J 


(53) 


where the quantity c is a function of X and v given by 


X— 1 

C (5 — 4r)X — (1 — 2r) 2 


(54) 


Appendix B contains the details of this calculation and others 
used in this section. It has been pointed out previously 
that N x and N„ must vanish. According to equations (53), 
this yields 


/ h V ijl 2ijl 

\2 — 2< y ~~cK~^R 


(55) 


Thus the strain rates e, and e 2 in the middle surface are 
related by the equation 


that is, 

z)>z 0 for positive K 1 

(51) 

z<z 0 for negative K j 


The surface z=z 0 separates the regions of loading and un- 
loading in the plate; a given part is in a state of loading or 
not according to which condition of equations (51) is satisfied. 

The criterion just developed must now be applied to the 
problem of buckling. As mentioned previously, the stress 
distribution of the buckled plate differs from the original 
state of pure compression by certain additional stresses 
c Tjjit , Vydt, and T zv dt. These new stresses are such that their 
total stress resultants must vanish and the moments pro- 
duced will be in equilibrium with the moment generated by 
the original compressive force in the budded plate. The 
vanishing of the stress resultants will lead to a formula for 
s 0 in terms of the constants of the material and the value of 
a 0 . Once this equation has been developed, a rather straight- 
forward computation will lead to the desired equation of 
equilibrium. 

The rates N x and N u of the stress resultants are defined as 


«!■ — ' — 2 


(56) 


From the definition of z 0 , equations (49), and this result, it is 
seen that 


Z 0 =(2-r) j+(2r-l) | 



(5 — 4r) 


(57) 


This result may then be substituted back into equation (56) 
to yield 


h V 2e 2 h 2z 0 h 
2~ Zo ) ~cK~ e(5- 4r) 


(58) 


Relation (58) can be solved for z 0 to yield the equation of the 
neutral surface which separates the regions of loading and 
unloading. Since h, c, and v depend only on the geometry 
of the plate and the applied compressive stress, z 0 will 
depend only on these quantities. It is more convenient to 
introduce a new' quantity defined by 

fo=x (59) 
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Then the quadratic equation for z 0 , equation (58), becomes a 
quadratic equation in to> namely, 


(Wo) 2 


4fo 

c( 5— -4k) 


= 0 


(60) 


There are two solutions to this equation, in general, but the 
only one which is physically realizable is 


to — fo +== «'b V<* 2 — 1 

where 

a=1 ~c( 5-4r) 


(61) 

(62) 


Equation (61) gives the desired formula for the neutral sur- 
face. 

When K<C 0, the procedure is exactly the same as that out- 
lined. Again the details are found in appendix B. The 
formula for the neutral surface is given in this case by 


fo=fo =— a— V<a 2 — 1 (63) 


Roughly speaking, the sign of K indicates whether the plate 
“buckles up” or “buckles down.” Consequently, the differ- 
ences in sign found by comparing equations (61) and (63) 
are quite natural. It would also be expected that the rates 
of change of the bending and twisting moments, as well as the 
resulting equilibrium equation, should be independent of 
the sign of K. This will be shown to be true. 

The rates of change of the bending and twisting moments 
can be computed now that z 0 or f 0 is known for i£> 0 and for 
0. The rates of change of the bending moments, M x 
and M v , are defined as 


Mx= 

P/ 2 , ^ 

<T x zdz 

J 

-ft/2 

1- 

% 

fft/2 

II 

SC 

*5 

dyzaz 

J 

-ft/2 J 


(64) 


be the same in both cases. The details are found in appendix 
B in which it is shown that 

M *=~ r2( f -C ) {KAl-cK2-vY]+K 2 lv-c5(2-v){2 V -l)]} 

(67) 

h 3 

My— — 12(1 p 2 ) {K l [v-c5{2-v){2 v -l)\+K t [l-cB(:2v-lY]} 

( 68 ) 

A '“i2(rb) [-1- &■] «*» 

The equation of equilibrium can be set up in terms of the 
bending and twisting moments and compressive load without 
reference to the stress-strain relations. This has already 
been done by Timoshenko (reference 1, p. 305) for the more 
general case of combined bending and tension or compression. 
His results may be applied to this special case of a simply 
compressed plate. With the present notation, the equation 
of equilibrium is 


b-M x b -M xy b 2 M y _ v b~w 

bx 2 “ bxby by- x bx 2 


(70) 


Timoshenko’s relation was originally written in terms of the 
actual bending moments and actual compressive stress re- 
sultant N x rather than the rates of the reduced quantities. 
Timoshenko’s equation can be differentiated with respect 
to time and divided by E 0 on both sides, so that equation (70) 
is the desired equation of equilibrium provided N x is defined 
as 

N z =(t<Ji/Eq (71) 

From equations (46), 

Ki = b 2 w/bx 2 
K 2 =b 2 w/by 2 


K l2 =2b 2 w/bxby 


where the moments are taken about the y- and ar-axes, 
respectively. The rate of change of the twisting moment 
M x v is given by 

. r h /2 

Mxu=— J_ h/2 T zv zdz (65) 


The calculation of the rates, equations (64) and (65), must 
be carried out separately for Kyi) and K<tb. It can be 
shown that the only quantity appearing in the final result 
which depends on the sign of K is the function 5 defined by 



[ l-ho + +l (fo + ) 3 ] for K>0 

[i+ir.-J (hr) 3 ] for k<o 


y (66) 


According to equation (63), = — fo + ; the numerical value of 

S obtained from equations (66) will therefore be the same in 
either case. Thus the expressions for M x , M y , and M zv will 


With these relations and equations (67), (68), and (69), 
equation (70) may be rewritten as 

E\\W x xxr~t' 2T-) \ 2 Vj xx yy -f- I-) 2 2'd']iyyn — (To hw x x (^2) 

where the subscripts denote partial differentiation with re- 
spect to the variable named and 


D n =D [1 — c5 (2-z/) 2 ] 

[1 — cS (2 — v) (2v- 1)] 
D22—D [1—cS (2v — l) 2 ] 
h 3 E 0 


D 


1 2 ( 1 — r 2 ) 


(73) 


The quantity D is the well-known flexural rigidity of the 
plate. Equation (72) resembles the equation for the buckling 
of an anisotropic plate (reference 1, p. 380). There is one 
important difference, however. In the case of an anisotropic 
plate, the coefficients D n , D u , and D 22 are constants of the 
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material; for the plastic case D n , D l2 , and D 2 2 are functions of 
o- 0 . In other words, the plate is anisotropic but this aniso- 
tropy is caused by and is a function of the compressive stress. 
Consequently, certain changes must be made in the standard 
procedure for calculating budding loads for anisotropic 
plates. Several examples illustrating this technique are 
given in the succeeding sections. Graphs of the quantities 


TYPICAL BOUNDARY CONDITIONS FOR THE FUNDAMENTAL EQUATION 
FOR A SIMPLY COMPRESSED PLATE 

The discussion of the boundary conditions for the buckling 
equation, equation (72), is facilitated by expressing the 
moment rates M z , My, and M xy in terms of the second deriva- 
tives of the deflection rate w and the stiffnesses D n , Du, and 
Z >22 introduced in equations (73). Thus, it follows from 
equations (67), (68), and (69) that 


Z?i/= 1 — c5(2 — v) 2 


Du=l — c5(2 — v) (2v — 1) 

and 

D 22 —\ — cd(2v — l) 2 


EqMx — DnWxx 1 D 1 2 (1 V )D\ Wyy 

E§M. V ~ \D \ 2 (1 v)D\w xx D 22 W VV *■ 

E a Mx y — ' D ( 1 v) V'zy > 


(74) 


as functions of X for i/=0.32 are given in figure 3 and the 
numerical values are listed in table I. 

TABLE I 


The following boundary conditions are typical in the budding 
of rectangular plates: 

(1) Simply supported edge at x— 0. The deflection rate 
w and the moment rate M x must vanish at this edge ; that is, 


VALUES OP Du’, Du', AND D r / 
FOR >' = 0.32 


X 

D,,' 

D\i' 

Dn' 

1.0 

1.0000 

1. 0000 

1.0000 

1. 2 

. 93030 

1.0149 

. 99680 

1. 4 

. 87447 

1. 0269 

. 99424 

1.6 

. 82847 

1. 0368 

. 99212 

1.8 

. 78973 

1.0451 

. 99034 

2.0 

. 75655 

1. 0522 

. 98882 

2.5 

. 69086 

1. 0662 

. 98580 

3.0 

.64178 

1. 0768 

.98355 

3.5 

. 60345 

1. 0850 

. 98179 

4.0 

. 57255 

1. 0916 

. 98037 

4.5 

. 54698 

1.0971 

.97920 

5.0 

. 52544 

1.1017 

. 97821 

G.O 

. 49104 

1. 1091 

. 97663 

7.0 

. 46459 

1. 1147 

. 97541 

8.0 

. 44365 

1. 1192 

. 97445 

9.0 

. 42652 

1. 1229 

. 97367 

10.0 

. 41226 

1. 1259 

. 97301 



Figure 3.— Graphs of Dn\ Du', and D 22 ' as functions of X for y=0.32. 


w= 0 

Du'Wxx \Di 2 (1 


—v)D] w„„=0 


I 


(75) 


for x=0. 

(2) Built-in edge at x=0. At this edge, the deflection 
rate w and the slope rate w x must vanish; that is, 


for x—0. 


w— 0 1 



(76) 


(3) Free edge at y— 0. For a free edge, the rates of the 
bending moment M u and of the equivalent shear load 
( — 2bMxu/<>x)-]-(c)My/dy) must vanish (reference 1, p. 300). 
Consequently, 


[D n —(l—v)D]Wxx+D 2 i Wyy =0 1 

\Di 2 ~\~ (1 D\tjjzxy DvnWyyu 0 


(77) 


for y~ 0. 

(4) Plane of symmetry at y= 0. If the buckled shape of 
the plate is symmetrical with respect to the plane y= 0, the 
rates of the slope w v and of the equivalent shear load 
(— 2DM X yl(ix) J r (DMy/dy) will vanish. Therefore, 

Wy = 0 

^ (78) 

[Di 2 + (1 — r) D] Wxxy + D 22 Wyyy = 0 J 
or 



for y= 0. Should these boundary conditions be given on 
other edges, the necessary changes in the formulas can be 
made easily. 

Several examples of the buckling of a simply compressed 
plate with various boundary conditions of the type just dis- 
cussed are considered in the next section. In all these ex- 
amples it is assumed that the plate is in the state of compres- 
sion previously described, a x — —o 0 , and that the edges 
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x—a and x=0, perpendicular to the direction of the compres- 
sive force, are simply supported. The other boundary con- 
ditions are specified for each example. Equations (75) are 
satisfied at x=0 and x=a if the deflection rate is written 
in the form 

w=j{y) sin ( mwx/a ) (80) 

where to is an integer. Thus a section of the plate in the 
buckled state, obtained by setting i/=Constant, is described 
by a series of sine waves; the integer to gives the number of 
half waves. Substitution of this expression for tb into the 
partial differential equation, equation (72), yields the ordi- 
nary differential equation 

( m-ir/a ) i D n f— 2 (rmr/a) 2 D 12 f n - f- D 22 f lv — <r 0 h (mir/a) 2 f (81) 


in which the Roman numerals denote the corresponding de- 
rivatives with respect to y. With 


P 


’=(¥) 

, <T 0 h / TO7r\ 

Dl2 \ & / 


2 -Pi 2 

Po, 

TOTtV P U P 2 2- 


P, 


(82) 


(83) 


P 12 P 22 

equations (82) can be written in the form 

/ lv - 2p 2 /‘ 1 +p 2 (p 2 — g 2 ) /= 0 

The general solution of this equation in f(y) is 

f(y)=c 1 cosh ry-\-c 2 sinh ry-\-cz cos sy+c 4 sin -sy (84) 

where c L , c 2 , c 3 , and c 4 are arbitrary constants which must be 
determined from the boundary conditions, and 


r= Vp(2+p) 
s=tJp(.Z-P) 


(85) 


possess solutions c u c 2 , c 3 , and c 4 which are not all zero yields 
an equation in r and s from which the critical stress <r 0 can 
be determined. 

SPECIFIC EXAMPLES 

Plastic buckling of a narrow strip; relation of present 
theory to beam theory. — It seems worth while to investigate 
how the present theory of the plastic buckling of plates is 
related to the Engesser-Yon Kdrman theory of the plastic 
budding of beams. Consider a rectangular plate of thick- 
ness h which is simply supported along x—Q and x—a, free 
along y=±b/2, and under the compressive stress a x — — <r 0 . 
As b~* 0 for a fixed value of a, the buckling condition for 
this plate might be expected to approach that of a beam 
which has the length a, is simply supported at its two ends, 
and possesses a rectangular cross section of height h and 
width b. 

The solution given by equation (80) automatically f ulfills 
the boundary conditions at x=0 and x—a. In addition, the 
symmetry condition along y= 0, equations (79), requires 
that the coefficients in equation (84) satisfy 


and 

Accordingly, 


rc 2 -\-sCi = 0 
r 3 c 2 — s 3 c 4 =0 

c 2 =Ci— 0 


( 86 ) 

(87) 

( 88 ) 


if the function f(y) is not to vanish identically. Boundary 
conditions, equations (77), for a free edge along y—b/2 fur- 
nish the equations 


[P 12 — (1— v)D] (c! cosh w+c 3 cos T0 + 


A 


: (r 2 Ci 


, rb „ sb 
cosh 2 ~ cos 


>< 


Equation (84) is the fundamental relation which must be 
studied for each particular case of budding from a state of 
simple compression. The boundary conditions of the type 
discussed lead to linear homogeneous equations for the con- 
stants c u c 2 , c 3 , and c 4 . The condition that these equations 


[Du+ (1 — v)D\ (rci sinh T j— sc 3 sin ^)+ 
D 22 (V 3 c x sinh y + s 3 c 3 sin ^ = 0 


0 


r (89) 


The boundary conditions corresponding to a free edge at y— — bj2 are fulfilled automatically if equations (89) are satisfied. 
Now these relations are linear homogeneous equations in Cj and c 3 ; if they are to have nonvanishing solutions, the determinant 
of the coefficients must be zero. Therefore, 

\D 22 r 2 - m ^f [D n - (1 -v)D) J cosh r | - ji D [D a - (1 -r)Z?]| cos | 

f{A^-5- 2 [A2+(l-r)0]}sinh^ 8 ^D 22 s 2 + m ^ [D n + (l-r)Z>] J sin ^ 

This determinantal equation can be reduced to the form 

(£ 2 —m 2 Tr 2 R 2 Q) 2 T) tan (T} 2j tm 2 ir 2 R 2 Qy^ tanh | = 0 


(91) 
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where 


R=b/a 

£ =rb=b^p(q+p) 
V —sb=bi/p(q_—p) 


(92) 


Q-[D u -(.l-v)D}IPnl 
If, for a fixed value of the span a, the width b of the plate 
approaches zero, R and hence £ and tj tend toward zero too. 
Accordingly, the functions £ tanh (£/2) and ij tan (jj/ 2) ap- 
pearing in equation (91) may be replaced by £ 2 /2 and i? 2 /2, 
respectively. For m=l, in particular, this relation can be 
reduced to 


kD»'-ED n '*-D a '*+<?D n '*=o 

in which the following notation has been used 

„2 


(93) 


k—12 


q-pa 


(1 r 2 ) 


-IV 2 ) /ZV 2 


(94) 


D n '=Du/D 
Diz =D 12 /D 

de=d 22 /d 

A— (Pn’Di i 

A detailed development of equation (93) is presented in ap- 
pendix C. If the plate under consideration buckles within 
the elastic range, D n ’ =D U ’ =D 22 ' =1] consequently, A=0 
and Q=v. According to equation (93) then, k=l — v 2 and 

ttWE 0 


o- 0 = 


12a 2 


(95) 


The compressive force P=a 0 bh under which the plate buckles 
is thus seen to equal 


ir 2 h 3 bE 0 -PEI 
12a 2 ~~~aF 


(96) 


where l — bh?j 12 is the cross-sectional moment of inertia. 
Equation (96) is Euler’s formula for the elastic buckling of a 
simply supported beam of span a. 


If the plate buckles after the compressive stress <r 0 has 
exceeded the limit of proportionality, the evaluation of equa- 
tion (93) becomes more difficult. Noting the relations given 
in equations (94), equation (93) can be transformed into 

k—D n '—Q 2 D 22 . (97) 

For r=l/2, this relation can be handled quite simply. A 
straightforward computation (see appendix C) yields the 
following result for the critical compressive stress <r cr —a 0 - 


7T 2 E 0 E 1 

3a 2 (1 + Vx) 2 


(98) 


The ci’itical buckling load P is given by 


P=a„bh 


t~E 0 I 

~hr 


4 

(i+Vx) 2 


(99) 


where again l=bh 3 /l 2 is the cross-sectional moment of inertia. 
Hence the reduced modulus E* is given by 


4-Ep = AE q E 

(i + Vx) 2 (V^+VA) 2 


( 100 ) 


according to the definition of X given by equation (13). This 
reduced modulus is identical with that found in the Engesser- 
Von K firman theory of the buckling of beams beyond the 
elastic limit. (See references 2 and 3.) Thus in the case of 
an incompressible material, v= 1/2, the critical load for a beam 
can be found as a limiting case from the theory of plates. 
It should be noted that, within the framework of beam 
theory, the critical stress is independent of the value of 
Poisson’s ratio v. 

This result is not necessarily true for materials which are 
not incompressible, v^lj2. The quantity 


k 12 cr 0 a 2 

1 — V 2 TT 2 h 2 E 0 


( 101 ) 


which is simply a constant multiple of the critical stress for 
a given plate, has been evaluated as a function of X for the 
case just mentioned, r=l/2, and for a material with r=0.32; 
the results are plotted in figure 4. Although the two func- 
tions agree at X=l, as previously proved, there is a marked 



Figure 4. — The quantity kftl—v*) as a function of X for two beams having the same shape but different values of v. 
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difference which, increases with increasing values of X. Con- 
sequently, two beams of the same shape but with different 
values of Poisson’s ratio will buckle at different critical 
stresses. 

At first sight, these results may seem somewhat contradic- 
tory since the analysis of plates developed herein is a generali- 
zation of the Engesser-Yon Karmln theory of the buckling 
of beams and yet this theory does not always appear as a 
limiting case of the present analysis. An explanation of this 
discrepancy is afforded by a closer examination of the posi- 
tion of the neutral surface z 0 as a function of the parameter X. 
The function z 0 depends on the plate thickness h and on v 
and X. The limiting process considered does not affect this 
relation and the neutral line in the beam is determined by 
the intersection of the neutral surface in the plate and the 
vertical plane y= Constant. The neutral line so determined 
need not coincide with the neutral line determined directly 
from the theory of beams. This difference is the basis of the 
discrepancies in k/(l — v 2 ) noted. 

According to equations (49), the function z 0 is given by 

(2 — v) e'i + (2v — 1) e ’2 

* n 

where 

K= ( 2 — v)K\-\- ( 2v — V)Ki 

The quantities e 4 and e 2 are the strain rates in the middle 
surface in the x- and ^-directions, respectively. For the 
small displacements studied herein, Ki is the rate of curva- 
ture of the line of intersection of the plane Constant and 
the middle surface. Similarly, K? is the rate of curvature 
of the line of intersection of the plane x= Constant and the 
middle surface. In general, the position of the neutral sur- 
face depends on the four quantities e h e 2 , K x> and iC and this 
holds in the transition from plate to beam. On the other 
hand, the Engesser-Von Karman theory of beams assumes 
that the position of the neutral line depends on e, and Ki 
alone. (See reference 1, p. 158.) The position of the neu- 
tral line will depend on the type of analysis used. Specifi- 
cally, it can be shown (see appendix C) that the position of 
the neutral line found by considering the intersection of 
the neutral surface and a plane y= Constant and the position 
determined by the Engesser-Von Karman method agree if 
and only if r=l/2. Therefore, it is to be expected that the 
transition from plate to beam is valid only under this con- 
dition. 

Buckling of a simply supported plate. — The critical load 
for a simply supported plate under edge thrusts, stressed 
beyond the elastic limit, has been discussed by Timoshenko 
(reference 1, p. 387) on a purely formal basis. The object 
of this section is the development of an analytical formula 
for the critical stress of such a plate on the basis of equation 
(72) and the subsequent comparison with Timoshenko’s 
results. 

According to equations (75), the boundary conditions for 
a plate simply supported at all four sides can be written in 
the form 

0 ) 

> ( 102 ) 

— D n w xx ~ [Z>i 2 — (1 — v)D]w yu —Q ) 


at a;—0 and x=a, and 
w=Q 

—D 22 w vv —[D 12 — {l—v)D]w xx —Q 


( 103 ) 


at y=0 and y=b. The boundary conditions at x=0 and 
x=a are fulfilled automatically if the function w is of the 
form given by equations (80) and (81). Equations (103), 
applied to the function f(y) at y=0, require that 

Ci+c 3 =0 

— £> 22 (r 2 c, — s 2 c 3 ) + [D 12 — (1 — v) D] (Cl + c 3 ) = 0 



Consequently, 


thus 


Ci—— c 3 
c 3 (r 2 -fs 2 ) = 0 


ci = c 3 = 0 


(105) 


(106) 


since r 2 -|-s 2 7^0 if buckling is to take place. 
The conditions at y=b imply that 

c 2 sinh rb+Ct sin sb= 0 
— D 22 (r 2 c 2 sinh rb—s 2 c t sin sb ) -f- 


y (107) 


Or)’^ 


- D (1 — v) ] (c 2 sinh rb + c 4 sin sb) — 0 


which are equivalent to 

c 2 sinh rb-\-Ci sin s6 = 0 
r 2 c 2 sinh rb—s 2 Ci sin sb = 0 


(108) 


Now these equations are linear, homogeneous equations in 
the two unknowns c 2 and c 4 ; they possess nontrivial solutions 
if and only if 


sinh rb sin sb 
r 2 sinh rb — s 2 sin sb 


= - — (r 2 -)-s 2 ) sin sb sinh rb — 0 


( 109 ) 


This condition will be satisfied if r 2 +s 2 =0, r— 0, or s=mr 
where n— 0, 1 , 2 ,.... The first possibility is untenable; 
the second implies, according to equations (108), that s=0, 
7T, 2 7r, .... If s=0, then from equations (108) r must also 
vanish and buckling will not take place. The buckling con- 
dition must t herefore be s—nw where n— 1, 2, 3, .... 

Since s= -y/piq—p) from equations (85), it follows from 
equations (82) that 


=p(q-p) 

rmr /Z?i 2 r l a oh (m^DuDv-Dn 2 mx lD 12 ~] ,, , 
~!TVD»l\D a \ a / D n D 22 aV AjJ (110) 
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With the notation of the previous sections for D n ', Dn, A/, 
and A, this relation can be written as 

n 2 ir 2 mir /D 12 ' r / <r 0 h ( mir \ 2 D 22 'A TOir /Z> 12 ' _ ( ,, , 

¥ ~ a V Dit L V DDE \a ) D n ' a V A/ J U ' 

The critical stress a cr can then be found by direct computation 

to be r - - ■ ■ 


&cr — <r 0 


t 2 DDE 

b 2 h 




. DA/ A 

mr~ 


(~y <“ 2 ) 


Now each of the terms in the first parentheses on the right- 
hand side is positive; consequently, the minimum critical 
stress will occur when n— 1. Thus the critical stress is given 

by 


/ 1 a /D 22 ’ 



DD 22 , A/mir 

' h 


m rV 

a / 


(113) 


In the elastic case, A= l and the functions appearing in the 
formula for the critical stress take on the values D n ,= 
D 1 2 , = D 2 2 , = 1, A=0; consequently, 


tt 2 D 

,<:r ~ 6 2 A, 



(114) 


This result is the same as that found for the elastic case 
directly (reference 1, p. 330). Timoshenko suggests (refer- 
ence 1, p. 387) the application of the following equation for 
the buckling of a simply supported plate compressed beyond 
the elastic limit, 

D (ct'W xr xz 2 ■yj Ot'lbxxyy ~\~Nz'W x x 0 (115) 

where here a=E*/E 0 , the ratio of the reduced modulus 
(defined in equation (100)) to Young's modulus. Under 
these assumptions, the critical compressive stress is found 
to be 

ir 2 D / la. b r \ 2 m i«\ 

b% \m b +m a ) (116) 

Although Timoshenko (reference 1, p. 387) gives only the 
critical stress for the half-wave number m— 1, his results can 
be easily extended to the form given in equation (116). 

It will be found more convenient to use a new quantity k 
rather than a 0T for comparison of the two theories. This 
parameter is a modified form of k defined in equation (94); 
more precisely, 


1 2Ccr¥ 

ir 2 h 2 E 0 


(1 — ***) 


t r c ,b 2 h 
ir 2 D 


(117) 


For the theory proposed in the present report, this quantity 
becomes 

* = - Dl 2 (mR V D?) +D ~ 2 ' Am2R2 ^ 118 ) 

where R=b/a, as before. In the elastic case, k takes the form 

j -(sb+» b )’ (u9) 

On the other hand, Timoshenko’s result can be written as 

I “(sB+”* ff ' fi ) ! (120) 

The differences in the two results appear more readily if equa- 
tions (118) and (120) are expanded in full. After some simpli- 
fications, equation (118) becomes 

k=D n 'm 2 R 2 +2D l2 ' + ^ 2 ( 121 ) 

while equation (120) takes the form 

k = am 2 R 2 -\-2-\/a J r^2'^2 ( 122 ) 


It is seen then that k is the same type of rational function of 
m 2 R 2 for both theories except for the fact that the coefficients 
are different functions of X. 

Unlike in the case of the cruciform section which is treated 
in the next section, the computations involved in evaluating 
lc as a function of X and mR are relatively simple since the 
coefficients D,,', D 12 ', D n ', and a can be tabulated once for all 
independently of the stress-strain law and the particular 
geometrical ratios under consideration. The critical values of 
k can be found with reasonable speed by the following pro- 
cedure. Curves of lc against X can be obtained for various 
values of the parameters m and R by evaluating equation 
(121) or (122), depending on which theory is used. These 
results, which form the bulk of the computation, do not 
depend on the stress-strain law but only on v. Conse- 
quently, the curves of k against X thus found are valid for all 
materials having the same Poisson’s ratio. On the other 
hand, k and X are related by a second equation which depends 
on the stress-strain law; namely, from equation (117), 


12 

k ~ ttWE 0 


(1-r 2 ) 


Since <r„ is a function of X, determined by the given compres- 
sive stress-strain law of the material in question, k also can 
be plotted as a function of X once the plate width-to-thickness 
ratio b/h has been fixed for a particular example. If these 
curves are plotted on the same sheet as those previously 
described, the intersections will give the critical value of 
lc, and hence the critical stress, corresponding to a given 
stress-strain law and given plate parameters b/h and R. It 
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should be noted that k depends on the ratio b/h in the plastic 
range, whereas this is not the case within the ordinary elastic 
theory. 

Since the quantity X enters into the computations for the 
second set of curves through the stress-strain relations, it is 
quite useful to represent the stress-strain curve analytically. 
An expression of the form 

e=^o+Aaf (123) 

has been fitted to the experimental data in such a way that 
the experimental and theoretical curves pass through the 
same initial and final points and the slopes at these points 
coincide. The fitted curve together with the experimental 
points is shown in figure 5. 

The computational program described has been carried 
out for the material given in figure 5 and the ratio of plate 
width to plate thickness fixed at b 2 /h 2 = 1000. The solid 
curves of figure 6 give k as a function of X with parameter R 
as defined by equation (121), each set of curves representing 
a different value of m. These curves can be applied to a 
plate with any ratio of b/h and any compressive stress-strain 
law provided r=0.32. The dashed lines, on the other hand, 
represent the relation between k and X given by equation 
(117); they correspond to b 2 /h 2 — 1000 and the stress-strain 



Figure 5— Stress-strain curve for e “ Jr Eo= 10,667,000, log A=— 42.873, and <o. =8.6127. 
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curve of figure 5. The points of intersection have been deter- 
mined and the critical values of k plotted in the usual fash- 
ion as functions of 1 jR in figure 7. Timoshenko’s results are 
also shown as well as those for the elastic theory. Since the 
critical stress is given by k multiplied by a constant, the 
interpretation of these curves of k against 1/R holds equally 
well for the critical stress as a function of 1/R. 

First, it should be noted that the plate will buckle in one or 
more half waves according to the magnitude of the ratio 
1 /R—a/b; this holds for plastic buckling as well as elastic 
buckling. Furthermore, it is necessary to consider only that 
part of the curve corresponding to a given value of m which 
lies below the intersections with curves belonging to adjacent 
values of to. Under these conditions, it can be readily seen 
that the critical stress obtained from the present theory lies 
between the results of the elastic analysis and Timoshenko’s 
formal procedure. More precisely, the elastic critical stress 
is higher than that predicted by the present theory, whereas 
Timoshenko’s buckling stress is lower. The transition from 
a buckling mode with a given number of half waves to the 
next higher occurs at practically the same values of 1/R in 
both plastic theories; whereas, this transition occurs at 
slightly larger values of 1/R in the elastic range. Finally, the 
locus of minimums of all three sets of curves is, to a high 
degree of approximation, a straight line. While this straight 
line is fixed for all values of b/h in the elastic case, it will 
shift in the plastic case as this parameter is changed. 
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oncluded. 

With a slightly different computational procedure, the 
results just obtained for the buckling of a simply supported 
plate by the present theory can be compared with Ilyushin’s 
solution of the same problem (reference 6). Ilyushin’s gen- 
eral equation for the buckling of a rectangular plate com- 
pressed in one direction is of the same form as equation (72) 
except for the difference in the coefficients. His fundamental 
relation can be written as 

k( [I* V) xxxx ~ k \3*' ti/xZI/y ~T" ‘22^/X)yyyy KoqWxX ( 1 24) 

where the coefficients D n *, D v >*, and D VI * are rather compli- 
cated functions of the stresses which will be discussed 
shortly. Since these results apply only to incompressible 
materials, the theory developed in the present report must 
be specialized to the case for which r=l/2. 

For the case of a simply compressed plate, the basic func- 
tions entering in the definition of the coefficients in equation 
(124) can be given in a slightly less complicated manner than 
that used by Ilyushin for the more general problem. Let 
E s denote the “secant modulus,” that is, quotient of the 
compressive stress divided by the compressive strain as 
obtained from a compression test for the material in question. 
Then Ilyushin’s function <u> can be shown to be 


co == 1 -■ 
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Figure 7.— -Comparison of values of k for values of 1/R with p= 0.32 as determined by different 

theories. 


where E 0 is Young’s modulus for the material. With Von 
K&rman’s modulus E* as given in equation (100), a new 
quantity k may be defined by setting 

k—E*/E 0 (126) 

The function i will he used to designate the following relation 
containing k and to, 


l/' = CO 




(127) 


The coefficients D n *, D 12 *, and D 22 * can then be shown to be 
of the form 


Z>„*=D(1-*) (l-i 1 ^^)' 

Z? 12 *=Z>( 1-*) 

D 22 *=D(l—\p) 


where D is the flexural rigidity of the plate defined in equa- 
tions (73). It should be noted here that Ilyushin’s coeffi- 
cients depend on the two moduli E and E„, whereas those 
appearing in the present theory depend only on the tangent 
modulus E or rather the ratio \=E/E 0 . Consequently, it is 
extremely difficult to carry out the larger part of the com- 
putations for k independently of the stress-strain law as can 
be done with the method presented herein. The numerical 
technique must therefore be changed somewhat. 

The boundary conditions for the simply supported plate in 
Ilyushin’s analysis are 

w= 0 

—D n *w xz — ^D 2 2*Wyy=Q 
at z=0 and x=a, and 



w= 0 

— 2 D 22 *w xx —D 22 *w ttt/ 

at y= 0 and y=b. These boundary conditions are com- 
pletely analogous to equations (102) and (103), provided 
Poisson’s ratio v is taken to be_l/2_in the equations (102) and 
(103). If the new quantities D u , D 12 , and D 22 are introduced 
in the following manner 



D n =Dn*/D 
D 12 =D i2 *!D *■ 

d 22 =d 22 */d _ 


(131) 


the critical parameter k can be found in precisely the same 
fashion as that used in developing equation (121). Con- 
sequently, 

k = DumH{ 2 + 2D l2 + (132) 

_ As has been previously noted, the coefficients D n , D 12 , and 
D 22 are functions of both the secant modulus and the tangent 
modulus with the result that specific reference to the stress- 
strain law in compression must be made in order to evaluate 
equation (132) readily. The critical quantity k has been 
computed as a function of a cr for the material shown in 
figure 5 with m and R appearing as parameters. For a 
given width-to-thickness ratio b/h, k can be determined by 
finding the intersections of these curves with the straight 
line defined by equation (117). This procedure has been 
applied to the case b 2 /h 2 — 1000 previously discussed. The 
results obtained from Ilyushin’s method and from the present 
method with r= 1/2 are shown in figure 8. Again consider- 
ing only those parts of the curves which lie below the inter- 
section points corresponding to consecutive values of m, 




PLASTIC BUCKLING OF A RECTANGULAR PLATE UNDER EDGE THRUSTS 


19 



Figure 8. — Comparison of values of k for values of 1/R with v=y> as determined by the present 
method and that of Ilyushin. 


certain general conclusions can be drawn. The minimum 
value of k as determined by Ilyushin’s method is practically 
constant for all values of m as has been seen in the other 
theories considered. On the other hand, Ilyushin's theory 
predicts smaller values of k and hence lower critical stresses 
than those of the present theory. Comparison with figure 7 
would indicate that Ilyushin’s critical stresses lie roughly 
between those of Timoshenko and the present theory. 
Finally, buckling occurs at practically the same wave number 
in both theories although the jump in form takes place at 
slightly smaller values of 1/R in the present theory. 

Torsional buckling of cruciform sections. — The torsional 
buckling of a cruciform section under compression can be 
studied by treating each flange as a simply compressed plate, 
simply supported at y = 0 and free at y=b. (See reference 1, 
p. 340.) As in the previous case, the deflection rate w will 
be given by equation (80) where /(y) is found from equation 
(84). The conditions for the simply supported edge (equa- 
tions (75)) at y — 0 become 


c 3 =— Ci ) 

S (133) 

/“( 0 )= 0 ) 

The last equation, implies r 2 Ci=s 2 c s . Since c 3 =-Ci and r 2 
and s 2 do not vanish, Ci=c 3 =0; consequently, 

f(y)=c 2 sinh n/+c 4 sin sy (134) 


Finally, the boundary conditions at the free edge y—b 
require, according to equations (77), that 

—(\-~cf) [£*12— (1 — r)Z>] (c 2 sinh ri+c 4 sin sb) + 

D 22 (r 2 c 2 sinh rb—s 2 c i sin sb)=0 

[Z?, 2 + (1 — v)D] ( rc 2 cosh rb + sc 4 cos sb) + 

I) 22 (r 3 c 2 cosh rb s ;i c 4 cos sb) =0 


-(“) 


(135) 


Now equations (135) are linear homogeneous equations 
in the two unknowns c 2 and c 4 . If these equations are to 
yield nonvanishing solutions, the determinant of the coeffi- 
cients must vanish ; that is, 


where 


a 1 
a 2 


61 ! 
h 


=0 


a l ={D 22 r 2 ~{in 2 ir 2 [D l2 —(l — v)D]/a 2 }) sinh rb 


(136) 


&i= — (Z? 22 s 2 +{m 2 ir 2 [Z7, 2 — (1 — v)D]/a 2 }) sin sb 


a 2 —r(D 2 2 r' 2 — {m, 2 ir 2 [D I2 -\- (1 —v)D]/a 2 }) cosh rb 
& 2 = — s(D 22 s 2 + {?« 2 ir 2 [Zli2+ (1 — v)D]/a 2 }) cos sb 


Except for a systematic interchange of the hyperbolic and 
trigonometric functions, the determinant in equation (136) 
is the same as that appearing in equation (90). The de- 
derminantal equation can then be reduced to 


s | Z? 22 r 2 — m J[- [Z>12 — (1 — v) D] | tanh rb — 

r {d 22 s 2 + T ~ [Z? I2 -(l-r)r>]J 2 tan sb (137) 

With the notation introduced in equations (92), this result 
can be rewritten in the form 

(£ 2 -mVB ! ® 2 — (,, 2 +mVtf 2 Q) 2 (138) 

Again £ and 17 can be determined as functions of the quantities 
D n ', D 12 ', D, 22 ', and k defined in equations (94). With some 
additional transformations, equation (138) can then be put 
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into a form suitable for solving for the quantity k as a 
function of X. Although the computations are somewhat 
tedious, the solutions obtained are completely independent 
of the stress-strain law. These results are represented by 
the solid lines of figure 9 and are marked for various values 
of the side ratio R. The dashed curves, on the other hand, 
represent the quantity k as a function of X as derived from 
equation (117). These curves have been computed for the 
stress-strain law shown in figure 5 and the values of the 
width-to-thickness ratio b/h as indicated in the graph. 
Given the two dimensionless quantities b/a and b/h, the 
critical stress is then determined by the corresponding 
solid curve and dashed curve. This point gives the desired 
value of k, from which the buckling stress is found by solving 
equation (117). The details of the procedure are found in 
appendix C. It should also be pointed out that the solid 
curves have been computed for the half-wave number m= 1, 
for it is shown in appendix C that the lowest value of k, and 
consequently the lowest critical stress, will be obtained for 
this value of the wave number. 



Certain general conclusions concerning the budding of such 
a section can be reached without reference to the stress- 
strain law by means of the solid curves of figure 9. These 
curves show that, for a given side ratio R, k is a decreasing 
function of X. This decrease is not the same for all values 
of the side ratio, however. For small values of R (large 
length- to-width ratios), k is almost constant. In fact, when 
R 5S 0.10, k is practically constant and the critical stress is 
the same as that found in the theory of elastic buckling of 
plates. The plastic effects are very pronounced for larger 
values of R and increase as R—> 1, that is, as the rectangular 
plate becomes more nearly square. 

Finally, the intersections of the solid and dashed curves of 
figure 9 yield k as a function of the side ratio R; these are 
replotted in figure 10. The top curve represents the elastic 
case (X= 1 ) and is independent of b/h. The other three 
curves correspond to the cases b 2 /h 2 = 300, 250, and 200, 
respectively. It should be noted that k for the elastic case 
is always greater than k for the plastic case; this difference 
is more marked for shorter (larger R) and thicker (smaller 
b 2 /h 2 ) plates. 



Figure 9.— Results of solving for A' as a function of X by using equations (117) and (138). 


Figure 10— Curves of k as a function of R from figure 9. 
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ENERGY METHOD 

In the preceding section rigorous solutions were obtained 
for several cases of plastic buckling of simply compressed 
plates. In cases for which the rigorous solution becomes too 
unwieldly or is not known at all, the energy ipethod provides 
a convenient means of finding approximate values of the 
budding loads. 

As far as the elastic buckling of plates is concerned, the 
energy method is well established. (See, for instance, 
reference 1, pp. 325 ff.) At first glance, it may appear 
somewhat doubtful whether this method can be extended to 
the plastic buckling of plates. That such an extension is 
legitimate, however, follows from the fact that equation 
(109) has the same form as the equation for the elastic 
budding of an orthotropic plate with flexural rigidities D n 
and Z >22 an -d torsional rigidity D n . Since the energy method 
is valid for this problem of elastic buckling and since the 
plastic buckling of a simply compressed plate is governed by 
the same equation, the application of the energy method to 
the plastic buckling of plates is legitimate. 

Under sufficiently small edge thrusts the flat form of the 
plate represents a stable equilibrium configuration. When 
the edge thrusts reach the critical value, however, this 
equilibrium becomes indifferent and the plate may assume a 
bent form. This transition from the flat to the bent form, 
that is, from one indifferent equilibrium configuration to 
another one, does not involve any energy input or output. 
Accordingly, the work done by the edge thrusts equals the 
flexural energy of the bent plate. The mathematical expres- 
sion of this principle is the basic energy equation: 

D//{[l-c5(2-.)W + 

2[v — c5(2 -v){2v — l)]w xx w vv - s r 

[1— c6(2r— 1) 2 ]w vv 2 + 2 (1 — v) w xy 2 }dxdy = trjbf f w x 2 dxdy 

(139) 

where the range of integration is the area of the plate; that is, 
O^x^a, O^y^b. As is easily seen from equations (46) and 
(74), the left-hand side of equation (139) equals 

—E 0 f f(M x K n —M xv K i2 +M v K 22 )dxdy (140) 

which represents twice the rate at which work is done by the 
bending and twisting moments. In the elastic range, c=0, 
this expression can be written in the familiar form 

Df f [ (w xx + w vv ) 2 — 2 ( 1 — v) (w zx w V y —w xv 2 ) ]dxdy (141) 

The right-hand side of equation (139) represents twice the 
ra'te at which work is done by the edge thrusts and has the 
same form as in the elastic case. 

The following expression for the critical load of the plate is 
obtained from equation (139): 

<t<a= 4 (i42 ) 


where 


A=dJ' {[1— cfi(2— v) 2 ]w xx -\-2[v— cS(2— v) {2v—Y)]w xx w yv - i r 
[1— c8(2v— l)]w V y 2J r 2(1 — v) w xv 2 }dxdy (143) 


and 


B-- 


-jj^dxdy 


(144) 


The right-hand side of equation (142) depends on the 
function w, that is, on the deflected shape of the plate. Now, 
any restriction which is imposed on the deflection rate w, over 
and above the boundary conditions discussed in the third 
section of the ANALYSIS, amounts to an increase in the 
stiffness of the plate and must, therefore, lead to an increase 
in the critical load. The critical load for a given plate is 
accordingly found as the smallest value which the right-hand 
side of equation (142) can assume for functions w possessing 
continuous partial derivatives of the second order and 
satisfying the boundary conditions for the deflection of the 
plate. 

In the clastic case, c=0, the right-hand side of equation 
(142) is independent of the buckling stress a 0 . Equation 
(142) therefore furnishes the buckling stress <r 0 which cor- 
responds to a given plate thickness h. In the plastic case, 
however, the numerator A of the right-hand side of equation 
(142) depends on the buckling stress a 0 through the quantities 
c and 5. Equation (142) should therefore be considered as 
an equation which determines the critical thickness h 
corresponding to a given buckling stress. If the expression 
for D, from equations (73), is substituted into equation (143) 
and the resulting form of equation (142) solved with respect 
to h, the following formula is obtained: 


/l2o- 0 (l — v 2 ) B 
V E 0 A' 


(145) 


where B is given by equation (144) and A' denotes the 
integral on the right-hand side of equation (143). 

The energy method mil now be applied to the buckling of a 
cruciform section previously considered. A suitable choice of 
the approximate deflection rate w is 


w=rj(S- 2?+?) (146) 

where £ =x/a and y—ylb. This function satisfies all the 
boundary conditions on the plate except at the free edge 
y—b. Nevertheless, this expression proves to be a satis- 
factory approximation as will be seen shortly. The necessary 
integrations indicated in the definitions of A and B (equations 
(143) and (144)) can be easily carried out in this case; for 
example, * 
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f J wjdxdy = ^ £ v 2 (a- d&r , = ^ (147) 

f f w xu 2 dxdy=~ R J\l-12^+8^+36^-48e+16^)d^ 


17 


"35 a 2 R 


(148) 


f f w x 2 dxdy=R £J\ 2 (1 -12£ 2 +8£ 3 +36£ 4 -48£ 5 -|- 16^)^ 

(149) 


17i? 
= 105 


where R=b/a. Thus 

B 17 a 2 R 2 

A' 168-S 2 [1 — c5 (2 — r) 2 ] -f- 102(1 — v) 


and 


17a 2 .R 2 

168tf 2 A/+ 102(1 -v) 
(150) 


f 168£ 2 A.' + 102(1 —v) 
k 17t r 2 


(151) 


where k is defined in equation (117). 



Fiquke II.— Curves of I obtained by the energy method and exact solution for >'=0.32 and 
various values of R. 


Since D n ' has been previously computed as a function of X 
(see table I), k may be readily evaluated from equation (151) 
as a function of X for various values of the parameter R. 
This has been done for the example r=0.32 considered in an 
earlier section. The results are shown in figure 1 1 in which 
the solid lines represent the solution obtained by the energy 
method and the dashed lines are those found from the exact 
solution. A brief inspection will show that the energy 
method, as applied here, gives a very good approximation in 
the technically interesting region of small values of R. The 
error increases as R becomes large and reaches a maximum of 
about 8 percent. It should also bo noted that for a given 
value of R the error is an increasing function of X. The 
energy method actually yields better results when the 
complete problem of determining the critical values of k 
associated with a given stress-strain law is carried out. 
These points are found from the intersections of the curves of 
figure 11 and the curves derived from the stress-strain law 
through equation (117) (fig. 9). Since the curves of figure 9 
arc monotonic, increasing functions of X with a slope angle of 
less than 90°, the error will be smaller than the original 
estimate. The results of applying this method to the cases 
treated in the previous section on exact analysis are shown in 
figure 12. The solid curves represent the solution by the 
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energy method, whereas the dashed curves are the exact 
solution. It is seen that the maximum error is of the order of 
2)4 percent of the value of k. This appears to be well within 
the order of accuracy of the theory itself. Finally, the 
approximate value of k is always greater than the exact 
solution as is usually found in the application of the energy 

method. 

Equation (151) also affords an easy means of studying the 
general character of the critical parameter k as a function of X 
for various values of R. It has already been pointed out that 
Du' is a decreasing function of X starting at Du' = 1 for 
X= 1 and approaching D u '=0.2413 for X— >°°. Thus for 
small values of R, the first term has little influence and k is 
approximately a horizontal line defined by the equation 


t 102(1 v) 
K ~ 17v 2 


(152) 


The elastic solution based on the energy method with the 
deflection given by equation (146) is 


r 168 jR 2 + 102(1 — v) 
17t r 


(153) 


and equations (152) and (153) agree quite closely for small 
values of R. Under these circumstances, then, the critical 
stress determined from the plastic theory will be the same as 
that obtained from the elastic formulas. 


Brown University 

Providence, R. I., March 26, 1947 

APPENDIX A 

DIFFICULTIES PRESENTED BY ILYUSHIN’S STRESS- 
STRAIN RELATIONS 

It has been pointed out in the first section of the ANAL- 
YSIS that certain objectionable discontinuities are present 
at the neutral surface when the state of stress and strain in 
a body is described by relations of the type used by Ilyushin 
(reference 6). More precisely, it has been shown that the 
equations for loading and unloading do not give the same 
increments in strain at the neutral surface; this difficulty is 
not present in the theory of plastic flow developed in the 
present report. Such problems always arise when an attempt 
is made to take unloading into account in a theory of plastic 
deformation. (See reference 10, p. 400.) 

The inconsistency is brought out more clearly if the whole 
problem is considered in terms of the state of stress and 
strain existing in the budded plate. Let 

fi=S2(o-() = r^— (Al) 

Then Ilyushin’s stress-strain relations for loading in the case 
under consideration (equations (30)) can be written in the 
form 


2E 0 e x = (1 + 12) (2cr z — <T„) 

2E 0 ty = ( 1 + 12) ( — a x + 2<r y) 

EoYiv = 3(l-|-fi)T IB 


(A2) 


Differentiated with respect to time, these become equations 
of the following type: 


2Eoi x — (1 -(-12) (2 i — <x„) 12 (2<r x <x„) 


(A3) 


In particular, for buckling from a state of simple compression, 
fT x 0*0 and (fy r xy — 0, 


0 = 


dti . 

dot ffi 


1 

2 


12' (2 <r x — <r„) 


(A4) 


where 12'=(/0/rfa- f . Accordingly, 

2E 0 e x = (1 "bli-t - vofiO (2 <Tz Gy) (A 5) 


The boundary between the regions of loading and unloading 
is a plane z— Constant. If the displacements are to be con- 
tinuous across this plane, i x , e v , and y xv , too, must be continu- 
ous. Moreover, according to Ilyushin’s theory the sign of 
— cr 0 e x serves as a criterion for loading and unloading. Since 
i x is continuous, it must vanish at the boundary between the 
regions of loading and unloading. Thus, 2 a x ~a v is con- 
tinuous, too, and 12 = 0 on the boundary. Denoting the 
values on the loading side of the boundary by the super- 
script + and the values on the other side by the superscript — , 
the strains at the neutral surface become 


Thus 


2E 0 i x = (1+12) (2<ti + — °V t 0 ' 

^ 2 & X (Jy 

= 0 

2E 0 i u =(l + Q)(-v x + +2« v + ) - 

= — oV + 2 <f !/ - 

•Eo7*!/ = 3(l+fi) T*» + 

=3t xv > 

ff x ~= (l+fi)<L + 

&y-=(l + tt)<ry + 

T X y = (1 “H fi) 



(A6) 


(A7) 


Since (1+0)>1, the rates of stress on the unloading side of 
the boundary are of the same sign as those on the loading. 
Furthermore, the absolute values of the rates of stress on 
the unloading side are greater than those on the loading side. 
While this prediction of Ilyushin’s stress-strain relations does 
not violate the equations of equilibrium, it does seem quite 
strange when compared with the usual notions of loading 
and unloading. The result is of a sufficiently startling 
character to call for direct experimental verification before 
a theory of structural stability in the plastic range is based 
on these stress-strain relations. 
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That the discontinuities implied in Ilyushin’s stress-strain 
relations are a source of concern to Russian scientists working 
in this field is seen from the following passage which is quoted 
from W. W. Sokolovsky’s recent paper (reference 14) in 
which the stability problem is specifically discussed: 

“In solving these problems, the Mises-Hencky theory was 
employed, under the assumption of the incompressibility of 
the plastic material. This assumption has a strong influence 
upon the value of displacements and does not permit satisfy- 
ing completely the conditions of the continuity of all components 
of stress and displacements on the boundary between the elastic 
and the plastic zones. An objective of recent work is to avoid 
these defects of the theory.” 

While the discontinuities appear to have been noticed, the 
source has not been traced correctly, for these discontinuities 
subsist even if the compressibility of the material is taken 
into account. They can be avoided only by replacing 
Ilyushin’s theory of plastic deformation by a theory of 
plastic flow of the type given in the present report. 

APPENDIX B 

DETAILS IN DEVELOPMENT OF FUNDAMENTAL PLATE 
EQUATION 

DETERMINATION OF NEUTRAL SURFACE 

The neutral surface can be determined from the condition 
that the additional stress resultants present in the buckled 
state must be in equilibrium. Rather than compute the 
rates of the stress resultants arising from ir x and <f„ separately, 
it is much simpler to consider the combinations ir x —vd v and 
h v —vir z . (A similar device is employed in reference 6, p. 
341.) For unloading, equations (40) may be used to yield 



The quantity c depends only on the properties of the material 
through v and on the compressive stress through A. It is a 
nondecreasing function of A, for 

dc (5 — 4r) A — (1 — 2v) 2 — (5 — 4r) (A — 1) 

d\ [(5 — 4r)A — (1 — 2r) 2 ] 2 

5— 4v — l + 4r— 4v 2 

[(5 — 4r)A — (1 — 2r) 2 ] 2 

4(1— v 2 ) 

[(5 — 4r)A — (1 — 2r) 2 ] 2 

^ 0 (B7) 


since — l^rAl/2. For A=1 (the smallest value of A), c=0; 
for A= co , c= 1/(5 — 4v). A similar expression for a y —va x can 
be found by the same technique; namely, 

<j v — vir x —i v — cK(z — 2o) (B8) 


The rates of change of the stress resultants N x and N v can 
be computed from equations (Bl), (B6), and (B8). Two 
cases arise depending on whether the region zf>z Q is a region 
of loading or whether s<s 0 is the region of loading; that is, 
whether iC>0 or whether K<f0. 

Case (1), Kf>0— Since iT> 0, the region zf>z 0 is the 
region of loading and z<fz a that of unloading. Consequently, 
equations (Bl) are valid for z^z 0 and equations (B6) and 
(B8) hold for 2 A s 0 . Now the rates of the stress resultants 
N x and N v are defined as 


phi 2 

N x ~ I a x dz 

./ -hi 2 


^h/2 
-hi 2 


( T ydfZ 


According to equations (Bl) and (B6) then 


(B9) 


On the other hand, equations (39) give for loading 


L is x ro-„+ (A l)vz 2 hy 


6 „=G v —va x 


A— 1 . , A— 1 . 


2 ~Gx~l ^ »v 

Transposition in the first of equations (B2) leads to 

• - A-1*,. . v 

V(Ty = t x 2 (2 Gx—Gy) 


(B2) 


(B3) 


which, according to equations (41), can be written as 
2(A— 1) 


v<Ty ~ ex (5 — 4f) A — (1 — 2v) 2 

With 

A— 1 


[(2— v)e x -\-(2v— l)e' v ] (B4) 


C— 


(5 — 4r)A — (1 — 2v) 2 


(B5) 


and K and z a as defined in equations (49), equation (B4) 
becomes 

h x — v<i tt =i x -\-2ciC(z— s 0 ) (B6) 


ph/2 . rh/ 2 

N x —vN v — e x dz+2cK (z—z 0 )dz (BIO) 

,/ —hj 2 J ZQ 

The strain rate e x =f—zK from equations (45); neither the 
function <=i nor K depends on z. Therefore 


and 


P /2 i x dz= P /2 dz-K( hl? zdz 

J-m J-h/2 J-ht 2 

=i 1 h (Bll) 

Nx-vNy^M-cK^-z^j (B12) 

In the same way, it can be shown that 

N y —vN x = f ,2 t u dz — cK f ,l/2 ( z—z 0 )dz 
J -A/2 J ?o 

=iji-\ck(^-z^ (B13) 

From equations (Bll) and (Bl3) it has been shown in the 

second section of the ANALYSIS that the position of the 

middle surface is given by that value of z 0 which is the solution 
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of equation (58). With Z 0 =2z 0 /h, equation (58) becomes 
equation (60); and with equation (62) the solution can be 
written as 

fo=a±V« 2 ~l (B14) 

Since c is a nondecreasing function of X with values in the 
range 0.5= 1/(5— 4r), the quantity a will lie in the range 

-oo^ajg-l (B15) 

Consequently, only the positive root in equation' (B 14) can 
be kept if the inequality |f 0 | = | 22 0 /fi,|^l is to be preserved. 
It will be convenient to denote this root by f 0 + ; that is, 

l'o + =fo=«+ V“ 2 — 1 (B16) 

This is the formula for the neutral surface given in the second 
section of the ANALYSIS. 

Case (2), K<^ 0, — The same technique, as used for case 
(1), with the exception of certain changes in sign can be 
employed in the case when 0. When 0, the region 
2 <\ 2 0 is the region of loading and z^>z 0 that of unloading. 
Consequently, equation (Bl) is now valid for z^ z 0 while 
equations (B6) and (B8) hold for Za- Thus N x —vN v 
becomes q 


<!+«•-; < B23 > 

and this can be transformed into 

? - ! + 2 [ I -i^E )] +1 “ 0 (B24) 

The quantity appearing in the brackets was defined as a in 
equations (64) ; thus the solution of equation (B24) is given by 

( B25 ) 

Again, the restriction that |f 0 | = l requires that only one of 
the signs in equation (B25) can be taken, namely, the minus 
sign. Thus the solution is that given by equation (63); 
namely, 

ro=-a-V^ II T=ro- (B26) 

where the symbol f 0 ~ has been introduced to distinguish this 
case from the root f 0 + previously found. Comparison of 
equations (B16) and (B26) shows that 

hf=-fo + (B27) 

In addition, f 0 + <0 for 2T> 0 and f<T>0 for K<C 0, as might 
be expected from the geometry of the situation. 


N 3 




A/2 

A/2 


e x dz+2cK 


if:;* 


dz (B17) 


which can be evaluated, as in the previous case, in the form 


N x —vN v =e\h—cK 

=iih —cK (j^+zo) (BIS) 

Similarly, 

Ny-vN x =e 2 h+\ K Q+ZoJ (B19) 

Again the rates of the reduced stress resultants N x and N v 
must vanish. This gives 



(B20) 


Equations (49) can be used to show that 




DETERMINATION OF RATES OF CHANGE OF BENDING AND 
TWISTING MOMENTS 

The moment rates can be computed now that z 0 or f 0 is 
known for each of the cases K^>0 and Since the rates 

of change of the reduced bending moments M x and M y are 
defined as in equation (64), the device of computing M x —vM y 
and M y —vM x may be used here. 

Case (1), A)> 0. — Application of equations (Bl) and (B6) 
yields 

M x —vM v = I 11/2 e x zdz+2cK ( m z{z~z 0 )dz 
J -A/2 J *0 

KJd , 2 0 3 h*Z 0 , 2o^ mom 

= W+ 2cK \2i ~^—\ s + TJ (B29) 

The quantity z 0 can be replaced by h£o + /2 — z 0 with the result 
M- vM v = [lY-clt (l -§ fo + +| fo +3 )] (B30) 

With 

25=l-|f„ + +|ro +3 (B31) 



(5 — 4r) 


and equations (B20) and (B21) yield 


(B21) 


the last equation becomes 

. . h 3 

M x —vM v =—j2 (K,-2 c8K) 


(B32) 


h . V 2 e 2 h 2zji 

2+V ~ ' cK~c(5- 4r) 


(B22) 


Equation (B22) is the analog of equation (58) which was 
found for the case K^> 0. It can also be solved by intro- 
ducing the nondimension al quantity defined as f 0 = 2 2 0 //L 
The quadratic equation for f 0 is 


Equation (B8) for a v —va x differs from equation (B6) for 
<r x —vdy only in that —K replaces 2 K and i v replaces i x in 
equation (B6). Without further computation, it can be 
seen that 

. . h 3 

My-vM x =-~ (Kt+c&K) (B33) 
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Equations (B32) and (B33) can be solved for the reduced 
bending moments with the result that 

- I2(1 ~V ) [Ki+vK 2 -{2-v)c8K] (B34) 

h 3 . . 

M v — — \2{^—v 2 ) [^+K 2 -{2v-l)c&K] (B35) 

The curvature K can be expressed in terms of K\ and K 2 ac- 
cording to equations (49); consequently, 

• h 3 

Mx= ~ 12(1~v 2 ) {Ki[l~cb{2~vY]-]-K i [v—c8{2—v) (2v— 1)]} 

(B36) 

~1 2(T— y 2 ) {KAv-c5(?-v) (2r — 1)]+-^2[1 — cS(2v l) 2 ]} 

(B37) 


APPENDIX C 


DETAILED COMPUTATIONS FOR SPECIFIC EXAMPLES 

BUCKLING OF A THIN STRIP 

According to the definitions of r and s, equations (85), and 
of p, equations (82), 

e-«>-2j) ! =2 (^J ^ (Cl) 

The terms of the second row of the determinant (equation 
(90)) can therefore bo written in the form 


r |a 2 s 2 +~t- [ 0 . 2 - (1 — r)Z?]| sinh ~ 
s [D 22 r 3 -^f [A 2 -(l-r)Z?]J sin ~ 


(C2) 


These relations are the same as equations (67) and (68) which 
were given in the second section of the ANALYSIS without 
proof. 

Case (2), K<^0. — This case can be handled in exactly the 
same manner as the previous one provided the integrals are 
split up as follows: 

M x — vM y — I 1 i x zdz+2cK\° z(z—z 0 )dz (B38) 
J ~k/2 J -ft/2 

= - tf + 2 c ^(t + S~y +£ ^) (B39) 

For K<C 0, z 0 —hg 0 ~/2 and 

M x -vM v =-^ l-\ f 0 - 3 +| fo“)] (B40) 

N ow fo _ = — f 0 + ; thus, 

(l-§ fo-+| fo")=(l+| fo +3 "! .(o + )=2a (B41) 

In other words, M x —vM y will be given by the same formula 
for AT<0 as for K'>0. Similarly M u — vM x will be given by 
the same expression in either case; and, consequently, 
M x and M v can be found from equations (B36) and (B37) 
independently of whether K<C _ 0 or K^>(). 

The rate of change of the reduced twisting moment M xy 
is defined in equation (65). According to equations (40) 
and (45) equation (65) may be rewritten as 


The determinantal equation, equation (90), is thus seen to be 
equivalent to 

« jz? 22 r 2 -^ [D 12 - (l-,)D] j' tan | 6 + 

r | D 22 s 2 + [D 12 ~ (1 - r) Z>] | 2 tanh = 0 (C3) 

With the definitions of equations (92), equation (C3) can 
be rewritten in the form of equation (91). The evaluation 
of the roots of the transcendental equation (91) can be simpli- 
fied further by introducing equations (94). A simple compu- 
tation will show that 



m 2 ir 2 R 2 , 

mVR 2 

D 2 


B -(V' 


IkDW 
j m 2 

- AZl 22 ,2 -(-Z?l 2 , ^ 

jkD 22 
\l m 2 

-AZV 2 — A/) 


(04) 


If, for a fixed value of the span a, the width b of the plate 
approaches zero, R and hence £ and n tend toward zero too. 
Accordingly, the functions £ tanh (£/2) and v tan (jj/ 2) ap- 
pearing in equations (C4) may be replaced by £ 2 /2 and rf/ 2, 
respectively. For m— 1, in particular, equation (91) can 
then be written as follows: 


( V kD 22 ' — AD 22 ,2 -\-D i2 — <2ZV) 2 (VHV- — AZV 2 — D l2 ) -f- 
(V*D«- — AD 22 ' 2 — D\ 2 QD 22 ') 2 (^JkD 22 ' — AD 22 2 -\-D l2 ) =0 


— 1 C ft/2 

M xv = - 2 (jqr^ J _ fc/2 ixyZ dz (B42) 

i r*/ 2 

= ~2(I+7)J_ A/2 (y-^ l2 )zdz (B43) 

Integration of the last equation furnishes 

^ xy== 12(l- v 2 ) ^ 2 ^ 12 (® 44 ) 

This result, too, holds independently of whether K ^> 0 or 

K< 0 . 


( 05 ) 

or in the form of equation (93) as 

kD 22 '—AD 22 ' 2 —D l2 ' 2 + Q*D 22 ' 2 =0 

Equation (99) for the buckling load when v= 1 /2 can be 
developed from the following considerations. According to 
equations (73) and (94), 


D n '=\-\c8 
D i2 =D 22 = 1 . 


(C6) 
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for r=l/2. From these results and equations (92), it is 
easily seen that Q= v, consequently, 


k = | (1 — 3cS) 


(07) 


when r=l/2. For this value of Poisson’s ratio, a straight- 
forward computation will yield 


X-l 


fo* 


3X 

X+l 

X— 1 

Vx — 1 
Vx+i 

X (X-(~ V3) 


(C8) 


Therefore, 


and 


cS- 


(i+Vx) 3 

(Vx-i)(Vx+3) 
3(l + Vx) 2 


k = ~(l-3c5)- 


4- (1 + VX) 2 

On the other hand, according to equations (94) 

9 <r CT a 2 


(C9) 

(CIO) 


k ~ir 2 h 2 E a 

where cr cr is the critical compressive stress. 

ir 2 E a h 2 1 

(T+VX? 

and the critical budding load P is given by 

tt 2 EqI 4 


(Cll) 
Consequently, 
(Cl 2) 


P = <r eT bh- 


(l+Vx) 2 


from equation (99). This is the desired result from which 
Von Karman’s equation follows. 

It has been pointed out in the third section of the ANALY- 
SIS that the neutral line of the beam and the neutral surface 
of the plate will coincide if and only if v=\/2. This can be 
seen quite simply from the following discussion. If v=l/2, 
then z 0 =t\/K according to equations (49). It follows from 
equations (C8) that 


fo + -= ; 


h 


Vx— i 


*v 1 


(Cl 3) 


The distance of the neutral surface from the lower surface is 
given by 


h . h 

2 +2o- Vx+T 


--h 


■y[E 


■y/Eo -f -yjE 


(Cl 4) 


This is the same result as that obtained in the direct analysis 
of the buckling of a beam. (See reference 1, p. 158.) 

Conversely, assume that the two positions of the neutral 
line agree; that is, 

Vx— 1 


fo“' 


Vx+i 


(Cl 5) 


where f 0 + is computed from equation (B16). Since equation 
(Cl 5) must be an identity in X, it must be valid for any 
particular value of X. For computational convenience, 
consider the case when X=4. Then 


Vx — 1 
Vx + 1 


19— 12r — 4v 2 


(Cl 6) 


3 (5 — 4v) 

where u— 8r 2 +12r— 23. Consequently, 
u 1 


Co -1 


3 (5 — 4v) '3(5 — 4r) 

_1 

3 


Vw 2 — 9(5 — 4v) 2 


and subsequent simplification leads to 
u— — 5(5 — 4v) 

Evaluation of u in terms of v yields 

{2v — 1) 2 = 0 

1 

v 2 


or 


(Cl 7) 

(C18) 

(C19) 

(C20) 


Therefore the two methods of determining the neutral line 
will agree if and only if r=l/2, and the transition from plate 
to beam will be valid only under this condition. 

BUCKLING OF A CRUCIFORM SECTION 

Transcendental equation (138) must be solved to determine 
the critical stress. A straightforward computation will 
show that 




D 12 

mVB 2 




mWR 2 / kD 


m‘ 


A 


2 0 


rCV^-^-^OJ 


(C21) 


Since 


>2 _2 

? V 


2mV 2 j? 2 Z> i: 

D 22 ' 


(C22) 
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equation (138) can be written entirely in terms of ij and tbe 
parameters B, Q,, D u ', and I) 2 /. For given values of to and 
R, tbe remaining parameters in the transcendental equation 
are functions of X alone. Thus n can be determined as a 
function of X by solving this equation by the usual iterative 
procedures. The quantity k can then be found as a function 
of X by solving the second of equations (C21). These results 
are independent of the stress-strain law of the plate material. 
On the other hand, 

k= iwE 0 (1 ~ vZ) (C23) 

that is, for a given material k can be written as a function 
of X once <x 0 is known as a function of X. 

It is a little more convenient to write the solution in terms 
of the quantity k rather than k. _ This parameter was defined 
in equation (117). Curves for k as a function of X for given 
values of to and R can be obtained in the following way. 
The quantity k can be computed, as previously described, 
by solving equations (138) and (C21); k can then be found 
from equation (117). These results are independent of the 
stress-strain law and consequently hold for any rectangular 
plate. On the other hand, k can also be computed as a 
function of X for a given stress-strain law according to 
equation (117). The resulting function will depend on the 
parameter i a //r. Once the side ratio R, the width-to- 
thickness ratio b/h, and the wave form to have been selected, 
the value of k corresponding to the buckling stress can be 
obtained by finding the intersection of the two corresponding- 
curves of k against X computed as just outlined. 

It has been pointed out in the third section of the ANALY- 
SIS that the lowest value of k, and hence the lowest critical 
stress, will be attained for to== 1. This can be seen easily 
from the following considerations. Figure 9 showed that 
for m=l and a fixed value of X, k increases for increasing R. 
With a simple change of variables, the solution for any value 
of m can be obtained from the solutions for m=l. Let 


R’ = mR k (C24) 



Then equations (138) and (C21), written in terms of k' and 
R', will be precisely the same as those for k and R when 


m— 1. Thus the curves of figure 9 can be used for any to 
provided R is replaced by R' and k by k'. For m>l then, 
k^>k' and R<CR' . In other words when m^l, the value of 
k for a given value of R and X will be larger than the cor- 
responding value of k for to — 1 . Since the dashed curves 
are monotonic, increasing functions of X, this implies that 
the critical value of k will be lowest for m=l. 
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